EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-EP/2003-020 
May 5, 2003 



Measurement of Exclusive p°p° Production 
in Two-Photon Collisions at High Q 2 at LEP 



The L3 Collaboration 



Abstract 

Exclusive p°p° production in two-photon collisions involving a single highly 
virtual photon is studied with data collected at LEP at centre-of-mass energies 
89 GeV < y/s < 209 GeV with a total integrated luminosity of 854.7 ptr 1 . The 
cross section of the process 77* — ► p°p° is determined as a function of the photon 
virtuality, Q 2 , and the two-photon centre-of-mass energy, W 71 , in the kinematic 
region: 1.2 GeV 2 < Q 2 < 30 GeV 2 and 1.1 GeV < W 77 < 3 GeV. 
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1 Introduction 



The exclusive production of p° meson pairs in two-photon collisions was studied by several 
experiments PQI2|- A prominent feature of the reaction 77 — > p°p° is the broad cross section 
enhancement observed near threshold, the origin of which is still not well understood [3] . Most 
experiments studied p°p° production by quasi-real photons, whereas only scarce data involv- 
ing highly off-shell virtual photons are available [2.. The interest in exclusive production of 
hadron pairs in two-photon interactions at high momentum transfer was recently renewed since 
methods for calculating the cross section of such processes were developed in the framework 
of perturbative QCD 0J. In these models, the exclusive process is factorisable into a per- 
turbative, calculable, short distance scattering 77* — > qq or 77* — > gg and non-perturbative 
matrix elements describing the transition of the two partons into hadron pairs, which are called 
generalized distribution amplitudes. 

This Letter presents results on the study of the two-photon reaction: 

e + e" -> e + e~77* -> e + e~p°p , (1) 

where one of the interacting virtual photons is quasi-real, 7, and the other one, 7*, is highly 
virtual. The squared four- momentum, Q 2 , of a virtual photon emitted by the incident beam 
electron 1 -' is related to the beam energy, Ef,, and to the energy and scattering angle of the 
outgoing electron, E s and 6 S by: 

Q 2 = 2E b E s (l - cos6 s ). (2) 

A scattered electron detected ("tagged") by the forward electromagnetic calorimeter used to 
measure the luminosity corresponds to an off-shell photon with a large Q 2 . The rate of such 
processes is considerably reduced as compared to production by quasi-real photons due to the 
sharp forward peaking of the angular distribution of the scattered electron. 

The data used in this study correspond to an integrated luminosity of 854.7 pb -1 and were 
collected by the L3 detector j^j at LEP. Of this sample, 148.7 pb -1 were collected at e + e~ 
centre-of-mass energies, y/s, around the Z resonance (Z-pole), with average \fs of 91 GeV and 
706.0 pb -1 at centre-of-mass energies in the range 161 GeV < y/s < 209 GeV (high energy), 
corresponding to an average of 195 GeV. This Letter presents the production cross section 
as a function of Q 2 in the restricted kinematical regions 

1.2 GeV 2 < Q 2 < 8.5 GeV 2 (Z - pole) (3) 

and 

8.8 GeV 2 < Q 2 < 30 GeV 2 (high energy), (4) 

and the two-photon mass interval 1.1 GeV < TV 77 < 3 GeV. The data are compared to Vector 
Dominance models jB] and to a recent QCD model [Zj. 

2 Event Selection 

2.1 Exclusive four-track events 

The reaction ((H), contributing to the process 

e + e~ — > e + e~7r + 7r~7r + 7r - , (5) 
Throughout this Letter, the term "electron" denotes both electrons and positrons. 
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is identified by a scattered electron and four charged pions measured in the L3 detector. Tagged 
two-photon events are accepted by several independent triggers: two charged-particle triggers 
jH] and an energy trigger demanding a large energy deposition in the luminosity monitor in 
coincidence with at least one track [§]. The combined trigger efficiency, as determined from the 
data itself, is (93.6 ± 1.3)% at the Z-pole and (97.9 ± 0.6)% at high energy. 

Single-tagged events are selected by requiring an electromagnetic cluster with energy greater 
then 80% of the beam energy reconstructed in the luminosity monitor, which covers the range 
25 mrad < 9 < 68 mrad of the electron scattering angle. At high energy the lower bound in- 
creases up to 31 mrad due to the installation of a mask to protect the detector from the beam 
halo. 

Event candidates are required to have exactly four tracks, with zero total charge and with 
a polar angle, 8, relative to the beam direction, such that | cos^l < 0.94. A track should come 
from the interaction vertex and have transverse momentum greater than 100 MeV. In addition, 
four-track events, incompatible with the pion mass hypothesis, are rejected using the energy 
loss information. 

Events containing muons are removed from the sample. A search for secondary vertices 
is performed, and events with reconstructed short-lived neutral kaons are rejected. An event 
candidate is allowed to contain no more than one electromagnetic cluster, with an energy below 
300 MeV and not exceeding 10% of the total energy of the four-pion system. 

To ensure that an exclusive final state is detected, the momenta of the tagged electron and 
the four-pion system should be well balanced in the plane transverse to the beam direction. 
Thus the total transverse momentum squared, p\ , including the scattered electron, is required 
to be less than 0.2 GeV 2 . This cut is also effectively a cut on the virtuality of the photon 
emitted by the untagged electron and thus ensures that the Q 2 variable, calculated from the 
measured parameters of the tagged electron using ((21), corresponds to the photon with the 
highest virtuality. 

2.2 Background estimation 

The contribution to the selected sample due to e + e~ annihilation is negligible. The background 
is mainly due to feed-down from tagged two-photon interactions producing a higher multiplicity 
final state, which is incompletely reconstructed. To estimate this effect two background-like 
data samples are selected. Firstly, we apply the same selection procedure discussed above re- 
leasing the charge-conservation requirement. Events of the types tt + 7t + 7t + tt~ and 7r + 7r~7r~7r~ 
are selected, in which at least two charged particles were undetected. Secondly, we select 
7r + 7T~7r + 7r~7r° events, requiring the 7r + 7r~7r + 7r - subsystem to pass the four-pion selection dis- 
cussed above without imposing the p 2 cut, and to contain in addition exactly two photons with 
effective mass in the range of ±15 MeV around the 7r° mass. We require p\ < 0.2 GeV 2 in order 
to select tagged exclusive 7r + 7r~7r + 7r~7r° events, and then consider only their 7r + 7r~7r + 7r~ sub- 
system to represent the background contribution. We assume that a combination of these two 
data samples gives a good description of the background from partially reconstructed events. 
Their p\ distributions, combined with the distribution of reconstructed Monte Carlo four-pion 
events, agree with the p\ distribution observed in the data. These distributions are shown in 
Figure Q for the restricted Q 2 -ranges and . 
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3 Data analysis 



3.1 Selected Sample 

In the region of four-pion mass W 77 > 1 GeV, 851 events are selected, 498 events at the Z- 
pole and 353 at high energy. The four-pion mass spectrum of the selected events is shown in 
Figure EK. 

The mass distributions of tt + 7t~ combinations, shown in Figures^ andEh, exhibit a clear p° 
signal, while the mass distribution of 7r ± 7r ± combinations, shown in FigureEli, has no resonance 
structure. In Figure Eb, the clustering of entries in the region of the crossing of the p° mass- 
bands gives evidence for a contribution of p°p° intermediate states. 

The p°p° production rate is determined as a function of Q 2 and W^. The resolution of the 
reconstructed variables Q 2 and W J7 is better than 3% and thus the event migration between 
adjacent bins is negligible. 



To estimate the number of p°p° events in the selected four-pion data sample, we consider 
non-interfering contributions from three processes: 



The data statistics is not sufficient to reach conclusions about contributions from subpro- 
cesses involving production of higher-mass resonances such as the /2(1270). Therefore in the 
present analysis we assume that the data is described by the processes © only. It was demon- 
strated that such a model provides a good description of exclusive four-pion production by 
quasi- real photons 0. 

Monte Carlo samples of the process (jBJ) are generated with the EGPC [E] program. About 
two million events of each process are produced for both the Z-pole and the high energy regions. 
The W 77 and Q 2 -dependence are those of the 77 luminosity function [TT] and only isotropic 
production and phase space decays are included. These events are processed in the same way 
as the data, introducing specific detector inefficiencies for the different data taking periods. 

For acceptance calculations, the Monte Carlo events are assigned a Q 2 -dependent weight, 
evaluated using the GVDM form-factor for both photons. Taking into account the detector 
acceptance and the efficiency of the selection procedure, the detection efficiency for each Q 2 
and W 71 bin is listed in Tables UH3J It is in the range of 10% — 25%, almost independent of the 
process. It slowly increases with Q 2 and slowly decreases with W 77 . 



In order to determine the differential p°p° production rate, a maximum likelihood fit to the data 
of the sum of the processes (jOJ) is performed in intervals of Q 2 and W^. The set, fl, of six two- 
pion masses, the four 7r + 7r _ combinations and the two 7r ± 7r ± combinations, provides a complete 
description of a four-pion event in our model of isotropic production and decay discussed above. 
This choice of kinematic variables allows to fully exploit the information specific to each one 
of the processes © and to obtain their contributions to the observed four-pion yield. For each 



3.2 Monte Carlo Modelling 



77 
77 
77 




iv tt it it , non — resonant. 



(6) 



3.3 Fit Method 
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data event, i, with measured variables Qi, we calculate the probabilities, Pj(f2j), that the event 
resulted from the production mechanism j. The likelihood function is defined as: 

3 3 

a = I1EW), E*i = i> ( 7 ) 

i j=l j=l 

where Xj is the fraction of the process j in the 7r + 7r~7r + 7r~ sample for a given Q 2 or W 77 bin 
and the product runs over all data events in that bin. The probabilities Pj are determined by 
the six-fold differential cross sections of the corresponding process, using Monte Carlo samples 
and a box method [13J. 

The fitting procedure is tested by applying it on various mixtures of Monte Carlo event 
samples from the processes (jEJ), treated as data. The contribution of the p°p° production 
process is always reproduced within statistical uncertainties, whereas, for small statistics test 
samples, large negative correlations, in the range of 60% - 75%, exist between the p°7r + 7r~ and 
7r + 7r~7r + 7r - (non-resonant) fractions. Both contributions are, however, necessary to fit the data. 
Therefore, in the following, only the p°p° content and the sum of the p°7r + 7r~ and 7r + 7r~7r + 7r~ 
(non-resonant) contributions are considered. 

To check the quality of the fit, the two-pion mass distributions of the data are compared 
with those of a mixture of Monte Carlo event samples from the processes (JHJ), in the proportion 
determined by the fit. The data and Monte Carlo distributions are in good agreement over the 
whole Q 2 and W 77 range; an example is shown in Figure El 

As pointed out in Reference El the tt + 7T~ system in the p°7r + 7r~ final state cannot have 
an isotropic angular distribution, since, in order to conserve C-parity, the angular momentum 
between the two pions has to be odd. We have verified that our results are insensitive to 
variations of the underlying angular distributions in the production model. In addition, a good 
agreement of the measured angular distributions of the data with those of the Monte Carlo is 
observed, as presented in Figure EJ 

4 Results 

4.1 Cross Sections 

The cross sections, Acr ee , of the process e + e~ — > e + e~p°p are measured as a function of Q 2 and 
Wyy and are listed in Tables together with the efficiencies and background contamination. 
The statistical uncertainties, listed in the Tables, follow from the fit. The differential cross 
section da ee / dQ 2 of the process (JTJ), derived from A<r ee , is also listed in TableHJ When evaluating 
the differential cross section, a correction, based on the Q 2 -dependence of the p°p° Monte Carlo 
sample, is applied, such as to assign the cross section value to the centre of the corresponding 
Q 2 -bin US]. 

To evaluate the cross section cx 77 of the process 77* — > p°p°, the integral of the transverse 
photon luminosity function, Ltt, is computed for each Q 2 and W 77 bin using the program 
GALUGA which performs exact QED calculations. The cross section <r 77 is derived from 
the measured cross section Acr ee using the relation Aa ee = L TT a^. Thus <r 77 represents an 
effective cross section containing contributions from both transverse (T) and longitudinal (L) 
photon polarizations: 

a 77 (jy 77 , Q 2 ) = a TT (W^, Q 2 ) + ea TL {W^ Q 2 ) , (8) 
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where &tt and otl are the cross sections for collision of transverse-transverse and transverse- 
longitudinal photons. The ratio of longitudinal to transverse polarization of the virtual photon, 
e, given, approximately, by the expression: 

2E S / E b 

e ~ i + (E s /E b y ' [J) 

is greater than 0.98, for our data. 

The cross section of the process 77* — ► p°p° as a function of W 7J , listed in Table |2] and El is 
plotted in Figure El together with the sum of the cross sections of the processes 77* — > p°7c + tt~ 
and 77* — ► 7T + 7r _ 7r + 7r~ (non-resonant). The statistical uncertainties of the sum of these two 
cross sections take into account their correlations. The p°p° cross section is dominated by a 
broad enhancement at threshold, already observed in the data at Q 2 ~ pQ and at moderate 
Q 2 (2] . The two cross sections are listed in Table d and plotted in Figure as a function of 

Q 2 ■ 



4.2 Systematics 

The uncertainty on this measurement is dominated by statistics. The uncertainty on the 
measured cross section due to the selection procedure, estimated by varying the cuts, is in the 
range 7% — 20%, affecting more the higher Q 2 region. Different form-factor expressions used for 
reweighting the Monte Carlo events and the variation of the acceptance contribute to an overall 
shift in the range 2% — 6%. The fitting procedure uncertainty mostly depends on the box size. 
It is estimated to be in the range 7% — 18% for the fits in Q 2 and in the range 8% — 30% for 
the fits in W 77 . 

To estimate the uncertainties of the background correction, the background determination 
procedure is performed using only the 7r ± 7r ± 7r ± 7r T or only the 7r + 7r~7r + 7r~7r° samples. A contri- 
bution in the range 6% — 11% is obtained. 

Collinear initial state radiation has little impact on the measurement since for 91% of the 
selected events the energy of the tagged electron exceeds 90% of the beam energy. 

All the contributions are added in quadrature to obtain the systematic uncertainties quoted 
in Tables CHSl 



4.3 Fits to the Data 

Figure 03 shows the result of a fit of the differential cross section da ec /dQ 2 to a form [T7j 
expected from QCD-based calculations [7j: 

daJdQ 2 ~ Qn{Q2+ i — >2)2 . (10) 

The fit is performed using the central value of the mass spectrum < W 77 >= 1.94 GeV. It 
provides a good description of the Q 2 - dependence of the data with an exponent n = 2.4 ± 0.3, 
to be compared with the expected value n = 2. Only statistical uncertainties are considered. A 
common fit of the data taken at the Z-pole and at high energy is justified by the almost constant 
values of the photon polarization parameter e, which determines the energy dependence of the 
cross section. 

In Figure |HK the data are fitted with two different form-factor parametrisations, leaving the 
normalization as a free parameter. A form suggested in Reference |H1 based on the generalized 
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vector dominance model (GVDM) [12], provides a good description of the Q 2 -dependence of 
the data whereas a steeper decrease is expected for a simple p-pole form-factor. 
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Table 1: Detection efficiencies, e, background fractions, Bg, and measured production cross sections as a function of Q 2 for 1.1 GeV < 
W 77 < 3 GeV for Z-pole and high energy data. The values of the differential cross sections are corrected to the centre of each bin. 
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Table 2: Detection efficiencies, e, background fractions, Bg, and measured production cross 
sections as a function of W 77 , for 1.2 GeV 2 < Q 2 < 8.5 GeV 2 , for the Z-pole data. The 
cross sections of the reaction e + e~ — > e + e~p°p°, 77* — > p°p° and of the sum of the processes 
77* — > p°7r + 7r~ and 77* — > 7r + 7r~7r + 7r~(non- resonant) are also given. The first uncertainties are 
statistical, the second systematic. 
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0.175 ±0.077 ±0.053 


0.31 ±0.14 ±0.09 
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0.012 ±0.022 ±0.009 


0.02 ±0.04 ±0.02 
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Table 3: Detection efficiency, e, background fractions, Bg, and measured production cross 
sections as a function of W 77 , for 8.8 GeV 2 < Q 2 < 30 GeV 2 , for the high energy data. The 
cross section of the reaction e + e~ — > e + e~p°p°, 77* — > p°p° and of the sum of the processes 
77* — > p n + n^ and 77* — > 7r + 7r~7r + 7r~ (non- resonant) are also given. The first uncertainties are 
statistical, the second systematic. 
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Figure 1: The pf distribution of the selected n + n^n + n^ data (points) in comparison with 
Monte Carlo distributions of four-pion events (open histogram) and the background estimated 
from the data (hatched histogram). The arrow indicates the selection cut on p\. 
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Figure 2: Effective mass distributions for the selected events: (a) Mass of the four-pion system, 
W 77 . (b) Correlation between the masses of two 7r + 7r~ pairs (two entries per event). The higher 
mass of each pair is plotted on the horizontal axis, (c) Mass of the it + it~ combinations (four 
entries per event), (d) Mass of the 7i ± n ± combinations (two entries per event). 
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Figure 3: Mass distributions of 7r + 7r~ combinations (four entries per event) for events with 
1.1 GeV < W 71 < 3 GeV in the fitted Q 2 intervals. The points represent the data, the hatched 
area the p°p° component, and the open area the sum of the p°7r + 7r~ and 7r + 7r~7r + 7r~ (non- 
resonant) components. The fraction of the different components are determined by the fit and 
the total normalisation is to the number of the events. 
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Figure 4: Comparison of data and Monte Carlo angular distributions: (a) | cos0 p |, the cosine 
of the polar angle of the p° with respect to the 77 axis in the 77 center-of-mass system; (b) 
I cos^ I, the cosine of the polar angle of the pion in its parent p° helicity-system; (c) A<f>, 
the angle between the decay planes of the two p° mesons in the 77 centre-of-mass system; 
(d) I cos# a b I, the cosine of the opening angle between the two n + directions of flight, each 
one defined in its parent p° rest-system. There are two entries per event in (a),(c) and (d) 
and four entries per event in (b). The points represent the data, the hatched area shows the 
p°p° component and the open area shows the sum of p°7r + 7r~ and 7r + 7r~7r + 7r~ (non- resonant) 
components. The fraction of the different components are determined by the fit and the total 
normalisation is to the number of the events. 
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Figure 5: Cross section of the process 77* — > p°p° and the sum of the cross sections of the 
processes 77* — > p°ir + iT~ and 77* — > 7r + 7r~7r + 7r~ (non-resonant) as functions of W 71: for (a) 
1.2 GeV 2 < Q 2 < 8.5 GeV 2 and (b) 8.8 GeV 2 < Q 2 < 30 GeV 2 . The points represent the data, 
the bars show the statistical uncertainties. The horizontal line for the highest W 71 bin indicates 
the upper limit of the 77* — > p°p° cross section at 95% CL: 0.26 nb for the Z-pole data and 
0.087 nb for the high energy data. 
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Figure 6: The p°p° production cross section as a function of Q 2 , for 1.1 GeV < W 77 < 3 GeV: 
(a) cross section of the process 77* — > p°p° and (b) differential cross section of the process 
e + e~ — > e + e~ p°p°. The points represent the data, the bars show the statistical uncertainties. 
The solid line in (a) represents the result of a fit based on the generalized vector dominance 
model (B| and the dotted line indicates the expectation for a p-pole form-factor. The line in (b) 
represents the result of a fit to a form expected from QCD calculations. 
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